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Wave Resistance : Some Cases of Three-dimensional Fluid 

Motion, 

By T. H. Havelock, F.E.S. 
(Received November 27, 1918.) 

1. Calculations of wave resistance, corresponding to a pressure system 
travelling over the surface, have hitherto been limited to two-dimensional 
fluid motion ; in those cases, the distribution of pressure on the surface is 
one-dimensional, and the regular waves produced have straight, parallel 
crests. The object of the following paper is to work out some cases when 
the surface pressure is two-dimensional and the wave pattern is like that 
produced by a ship. A certain pressure system symmetrical about a point is 
first examined, and more general distributions are obtained by superposition. 
By combining two simple systems of equal magnitude, one in rear of the 
other, we obtain results which show interesting interference effects. In 
similar calculations with line pressure systems, at certain speeds the waves 
due to one system cancel out those due to the other, and the wave resistance 
is zero ; the corresponding ideal form of ship has been called a wave- free 
pontoon. Such cases of perfect interference do not occur in three-dimensional 
problems ; the graph showing the variation of wave resistance with velocity 
has the humps and hollows which are characteristic of the resistance curves 
of ship models. 

Although the main object is to show how to calculate the wave resistance 
for assigned surface pressures of considerable generality, it is of interest to 
interpret some of the results in terms of a certain related problem. With 
certain limitations, the waves produced by a travelling surface pressure are 
such as would be caused by a submerged body of suitable form. The expres- 
sion for the wave resistance of a submerged sphere, given in a previous 
paper, is confirmed by the following analysis. It is also shown how to extend 
the method to a submerged body whose form is derived from stream lines 
obtained by combining sources and sinks with a uniform stream; in par- 
ticular, an expression is given for the wave resistance of a prolate spheroid 
moving in the direction of its axis. 

2. Take axes 0^ and Oy in the undisturbed horizontal surface of deep 
water, and 0:2; vertically upwards, and let ^ be the surface elevation. For an 
initial impulse symmetrical about the origin, that is for initial data 
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where m^ = x^ + y^, the velocity potential of the subsequent fluid motion, and 
the surface elevation are given by* 

1 f^ 
(j> = -\ f (/€)e''^ Jo (fcm) cos (KYt) /cd/c J (1) 

f = -~ -i P/ (k) Jo (fc^) sin (fcYt) K^dfc, (2) 

9PJ0 

where V^ = gf/Cj and 

/(k) = [^F (ot) Jo (koc) u da. (3) 

We have assumed that it is permissible to use the integral theorem 

poo poo 

F (ny) = I Jo (fcm) /cd/c] F (a) Jo (/cot) a doc. (4) 

Jo Jo 

We obtain the effect of a pressure system moving over the surface with 
uniform velocity e in the direction Ox by integrating (1) and (2) after 
suitable modifications. We replace t hj t'-r and x by x--ct, and integrate 
with respect to r over the time during which the system has been in motion. 
We shall limit the present discussion to the case when the system has been 
in motion for a long time, so that if we take an origin moving with the 
system a steady relative condition has been attained. In this case, with a 
moving origin 0, we replace x by x + c^o and ^ by -z^ in (1) and (2), and obtain 
the required resultsf 

/•OO 1*00 

pcl> = e-^'^'^duX f{K)e'''^JQ[K^y{{x^\■cuf-\-y^}]Qo^{KYu)KdK, (5) 



fOO TOO 

e-i'^^'dici f(fc)Jo[fc^{(x + cuy + y^]Bm(fcYu)fcWdK, (6) 

where /(/c) is obtained from the assigned pressure distribution p = F (m) by 
means of (3). 

The introduction of the factor exjo(—fjLu/2) is familiar in these problems 
and needs only a brief explanation. It may be regarded as an artifice to 
keep the integrals determinate, it being understood that ultimately fju is to be 
made infinitesimal. Or, again, it ensures that the solution is the fluid motion 
which would establish itself eventually under the action of dissipative forces, 
however small. 

In the steady motion with which we are concerned, we may imagine a 
rigid cover fitting the water surface everywhere and moving forward with 
uniform velocity c. The assigned pressure system F(t3r) is applied to the 

* Lamb, * Hydrodynamics,' p. 425. 

t 'Roy. Soc. Proc.,' A, vol. 81, p. 417 (1908). 
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water surface by means of this cover; hence the corresponding wave resist- 
ance is simply the total resolved pressure in the direction 0^.* With the 
nsual limitation that the slope of the surface is everywhere small, this leads 
to 

E = |F(,.)^iS, (7) 

taken over the whole surface. 

The evaluation of the steady wave resistance for an assigned pressure 
F (m) is to be carried out by means of (3), (5), (6), and (7). However, we 
may obtain simpler expressions before applying them to particular cases. 

3. For this purpose we analyse the wave disturbance (5) into simple con- 
stituents, in fact into one-dimensional disturbances ranged at all possible 
angles round Ox, the line of advance. We have 

7rJo[/c{(«4-c?/)2 4.2/2j.i/2] -. g^*'^(^+<^w)cos<^cos(/c3/sin^)^<^. (8) 

Jo 

Substitute in (5) and we can now carry out the integration with respect to u ; 
for we have 







^-^fiu ^iKCU COS 4. gin (fcYu) du 

= («ccos(5!) + /cY + |-/i'2l)""i"-(/i:ccoS(^— /cV"f |/tii)~^. (9) 

We simplify this expression further by using the fact that /x was 
introduced only to keep the integrals determinate, and is eventually to be 
made infinitesimal ; we can therefore reject terms in /ju which are super- 
fluous for this purpose. The process receives its justification in the course of 
the analysis. This being understood, we can use, instead of the right-hand 
side of (9), the expression 

— 2 (V/c^)sec^^/{/c— /cosec^(j6-i-'i(/tt/c)sec^}, (10) 

where /cq = g/c^. Using these results in (6), and making a slight trans- 
formation, we can express the surface elevation in the form 

? = n S®^ ¥># '^/('^) i 3 , ... / . T 

g—tKix COS <l)+y sin <l>) ~\ 

In (11) we have the surface elevation analysed into plane wave con- 
stituents, each element moving in a line making an angle <p with 0^. Carry- 
ing out the integration with respect to /c, we can express each constituent 

* * Eoy. Soc. Proc.,' A, vol. 93, p. 244 (1917). 
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in terms of a simple harmonic wave in rear of the line x cos <^ + 2/ sin ^ = 0, 
together with a* disturbance symmetrical with respect to this line. We 
might continue the discussion for general types of pressure distribution, 
provided the functions are such that certain transformations may be used ; 
however, it is simpler to study in detail a few cases for which the conditions 
are all satisfied. 

4 At this stage it is convenient to specify the system 

p = ^{^)==Afl{f' + r,''r^ (12) 

where A and /are constants. Using (3) for this case, we have/(/c) = Ae""-^. 
Eeturning to the expression (11) for the elevation, we consider the element 
making an angle (f> with Ox. We change to axes Ox\ Oy' given by 
x' ■=^ X cos^-f-^/ sin <^, ^' = y cos i^—x sin^<3f>. The integral with respect to k 
then becomes 

A I tce-''^ \ - --{• - -. \dK. (13) 

Jo L^-— /^osec^^ + '?'(//'/ c) sec ^5) /c—fCoseG^^'—i{/ui/G)seG^j 

This integral can be transformed by contour integration ; as it is of a 
type familiar in plane wave problems^ we write down the results when these 
have been simplified by making /jl zero after the transformation has been 
carried out. We have for the value of (13), 
4:7r/coA sec^ ^e-''o/sec2</> ^{j^ (^^^^' ^^^2 ^^ 

. ^. r°°A:osec^(f)Cos/m + msin/'??^ ^. , „ , ^ 
Jo m2 + /co^sec^<^ 



2A 



00 



.osec^./,^cos/m + msin> ^_„.,^^^_ f^^ ^, ^ ^ ^^^^ 



nr + K(f sec* (^ 

From (11) and (14) we could now write down the elevation fas the sum 
of the constituents for all values of ^ in the range from — 7r/2 to 7r/2. The 
first term in (14) represents simple waves in the rear of the corresponding 
wave front x cos ^-{-y sin <^ = 0; hence the integration of this term 
would only extend over elements for which the assigned point {x, y) was in 
the rear of the wave front. The other terms in (14) represent a disturbance 
symmetrical with respect to the wave front, and diminishing with increasing 
distance from it. We shall not write down the expressions, as we do not 
intend to examine the wave pattern in detail. From the definition in (7), 
it follows that we can evaluate the wave resistance E by considering first a 
simple constituent of the elevation and then summing with respect to </>. 

Since the pressure system is symmetrical with respect to the origin, the 
symmetrical local disturbance in (14) gives no resultant contribution to R ; 

* Compare, for example, *Roy. Soc. Proc.,' A, vol. 93, p. 524 (1917). 
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also the part due to the regular waves in rear of the wave front 
X cos <j) + y sin <^ = is given by 

. o A Q I ^0.2^ P 7 / f^ A/eos(/^o^'sec2(f)) , , 

Collecting these results we have, from (11) and (15), 

K = (4:7rfgp) AW f sec^ ^^-^'^o/sec^ ^^^ (16) 

Jo 

We may express E in terms of known functions in two convenient forms. 

-A. 

If Wit, ^ (a) is the confluent hypergeometric function defined, under certain 
conditions, by^^ 

Wfo.(«)=: j;^|^^£M-*+»-i(l + ^/a)*+'»-ie-«6?t*, (17) 

and if K^ (a) is the Bessel function for whichf 

fco * 

g-acosh »/. (.Qgt^ r^^ ^^^ (18) 

we find, after some reduction, that 

E = {nr^^^l^gpp) A^a^/^^--/^ Wi, 1 (a) (19) 

= (7r/%p/')A'^'^~"^MKo(^/2)-(l + l/a)Ki(a/2)}, (20) 

where a = 2/<:o/ = 2^//c^. 

In a previous paper,J the same function of velocity, except for the con- 
stant factors, was found for the wave resistance of a submerged sphere ; the 
result was given in the form (19), and a graph was drawn to show E as a 
function of c. The resistance rises to a maximum in the neighbourhood of 
c = x/{gf)j and then falls asymptotically to zero. 

Although there are few tables available for the functions Kn in general, Ko 
and Ki are given in ' Funktionentafeln ' (Jahnke u. Emde) under the form of 
{i7r/2) Ho<^^ (ix) and (7r/2) Hi<^) (ix) respectively. 

5. Eeference has been made to the wave resistance of a sphere submerged 
at a depth/, large compared with the radius a ; this was calculated directly 
as the resultant horizontal pressure on the sphere. The connection with the 
present analysis is easily shown. 

In the paper referred to, the approximate solution for a submerged body 
was found directly, following Prof. Lamb's method for a cylinder. It is con- 

* Whittaker and Watson, * Modern Analysis,' p. 334. 
t Grey and Mathews, * Bessel Functions,' p. 90. 
t *Roy. Soc. Proc.,' A, vol. 93, p. 530 (1917). 
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venient to repeat here the expressions for the velocity potential and surface 
elevation due to a cylinder and to a sphere, putting them into the same 
notation for purposes of comparison. We have, for a cylinder 



«/> 



ca^x ca?x 



«' + («+/)' «' + («+/)' 



'{■2ca?\ 6~'^"/2 ^^ g-K if-z) gjjj k{x-\' m) sin {kYu) KSfdic, (21) 

Jo jo 



f = 2 a2//(a;2 4-/2) - 2 aM ^-^^/^ ^^^ ^-^z ^os «: (aj + c%) sin {kYv) kY clfc, (22) 
and for a sphere 



<#> = 



c<x% ccA 



r^ r 





•00 



"e--</~^) Jo {k^ {{x-\-cuf-\'y^}] sin (/cVi^) /cV^/c, (23) 



f = a^fjix^-^-y^^ff^ 



a? 



"^~mt2 ^^^ [ "e~«/ Jo [a:^/ {(^ + ci*)2 + y^)1 sin (a;Vw) /c^V^/c. (24) 

Jo 

These expressions satisfy the conditions at the free surface, namely, 
cd(p/dx-\'g^ = and d^/dz = ed^/dx, when fju is made infinitesimal. Oppor- 
tunity has been taken to correct an obvious mistake in sign in the expres- 
sions for the sphere ; in the former paper, the last terms in (23) and (24) 
were given as positive instead of negative. 

Eeturning to the comparison with § 4, consider the expression (24) for the 
surface elevation due to a submerged sphere. The first part represents a 
disturbance symmetrical about the origin, due to a doublet at the centre of 
the sphere, together with an equal opposite doublet at a point a height/ 
above the free surface. Compare now the second term in (24) with the 
surface elevation given by (6) when the pressure system is (12), so that 
/(/c) = Ae'^f The two expressions are identical, with a suitable relation 
between the constants ; we must have a^ = A/gp, or the corresponding 
moment of the doublet is Ac/2gp. We have then two related problems. 
For the submerged sphere the pressure is constant at the free surface, and 
the surface elevation consists of the two parts in (24) ; the wave resistance 
depends upon the supply of energy needed to maintain the waves contained 
in the second part of (24), and this energy is supplied through ^the work of 
the pressure at the surface of the sphere. On the other hand, for the travel- 
ling surface pressure, 

p=gpaY/{f' + m'f^ (25) 
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the surface elevation is the same as the second part of (24) ; hence the same 
supply of energy is needed, and is obtained in this case from the work of the 
applied pressure. Thus we may assume that the wave resistance is the same 
in the two problems. 

From (19) and (25) we have 

E=i7r%p/'-%«a3%-«/3Wi,i(^X (26) 

which agrees with the value for the sphere given in the previous paper. The 
connection between the wave patterns of a submerged body and of a certain 
surface pressure has been pointed out by Dr. G. Green in a recent paper,* in 
which the correspondence is developed from a different point of view. In the 
following analysis we deal only with combinations of simple pressure systems 
(12), and the corresponding submerged body can be found from a similar 
combination of doublets, as in the preceding case. 

6. The foregoing results can be generalised for other symmetrical forms of 
local pressure distribution, provided transformations such as are used in (4) 
and (14) are applicable. Assuming this, it appears, from the analysis of § 3 
and § 4, that for a pressure system p = F (^sy) we have 

E = (4 irlgp) ko^ ^ 'sec^ ^ {/(% sec^ ^) }2 J^, (27) 

Jo 

where /(/c) is given by (3). 

7. Some points of interest in the theory of wave resistance can be 
illustrated by combinations of the simple type (12). We consider first two 
equal systems, at a distance 2 A apart, and advancing in the direction of the 
line of centres ; that is, 

p = A//(/2 + m'fl' + A//(/2 + m^f/\ (28) 

where mi^ = (a? — hf + y^ and ^2^ ^ {x + hf + f. 

Writing, for the moment, jpi and p2 for the two component systems, and 
?i, ?2 for the surface elevations which would be caused by these systems 
acting separately, the waves due to the combination are given by ^1 + ^2, since 
we neglect, as usual, the squares of the fluid velocities. It follows from the 
definition in (7) that the wave resistance is the sum of four parts, Ri, E2, E12, 
and E21. Here Ei is the resistance due to the pressure pi acting on the waves 
produced by pu E12 is that due to pi operating on the waves caused by p2f and 
similarly for R2 and E21. 

It follows from § 4, that 

El = E2 = (47r/^p) AW r''sec^(^6""W«««^*#, (29) 

* a. Green, * Phil. Mag./ toI. 36, p. 48 (1918). 
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The terms E12 and E21 represent the interference effects. Let Bi and B2 
be the centres of the two systems pi and J02. To calculate E21, consider a 
constituent plane wave-front through Bi ; take this line as a new axis O'y' , 
and a perpendicular line through B2 for the axis 0'x\ as in %. 1. 




Then, corresponding to the expression (15) in § 4, we have as an element of 
E21 the quantity 

4^«o^Asec3,^e-/-»4°°%'r ff ,7 ^"If/ ff ^^.^ (30) 

The similar element in the value of E12 is the expression (30) when we 
have replaced ^' + 2Acos^ by .-z;' — 2Acos^. Adding the two elements and 
carrying out the integration with respect to y\ we have, as an element of 
E124-E21, 

'^ ^ J_oo(aJ +2Acos^)2-|-/^ 

= 8 7rWA2 sec^ ^^-SKo/sec^^ cos (2fCoh sec <^). (31) 

Eeplacing from (11) the proper factor and integrating with respect to <}>, 
we have 

E12 + E12 = (STr/gp) Ko^A^ f sec^ ^^-SKo/sec^^ cos (2fCoh sec <^) ^0. (32) 

Finally, from (29) and (32), the total wave resistance E is given by 

E =: (16 tt/^p) aw r sec^ ^^-s^co/sec^ cos2 (fcoh sec ^) d(f>, (33) 

We can express E in series of known functions by expanding cos(2Aro>^ sec ^) 
either in powers of fCoh, or in Bessel functions Jn(2/CQh); however, as these 
series involve either 'W]c,m{^/Cof) or K„(a:o/), they are of no use for numerical 
calculations. 

It is not difficult to calculate numerical values directly from the integral 
(33) for given values of the constants. To obtain a graph showing the 
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variation of the resistance E with the velocity c, the following method is 
sufficiently accurate, at least for illustrating the main features. 

Take as a definite example, h = 2/; then, writing a for 2^/c^, we require 

•7r/2 

sec^ <^e~«sec2<^ (,Qg2 ^^ gg(, 0^ ^0 fQj. various values 







to calculate the value of o? 

of a. 

The integrand can be obtained without much trouble, and it was found 
sufficient to calculate its value at intervals of 10° throughout the range from 
to 7r/2 ; the mean value was found from half the sum of the initial and 
final values together with the sum of the intermediate ones. In the course of 
these calculations, we have material for obtaining the value of 

r7r/2 
. 

by the same method ; but this integral is equal to 

la^^-«/2{Ko(a/2)~(l + l/^)Ki(a/2)}, 

and we can find its value also from the tables of Kq and Ki mentioned in | 4. 
By comparing results we obtain some idea of the accuracy of this method of 
numerical integration. The calculations can be lightened for present 
purposes by choosing, from general principles, values of a which correspond 
to important points on the graph. 

By this method we obtain values of E for different values of c, for this 
particular case. The result is shown in the full curve in fig. 2 ; the scale for 
E is arbitrary, the unit for c is the velocity \/ {gf)- The dotted curve is a 
mean curve, and is equal to Ei-i-E2 in the notation of this section ; that is, 
it represents the sum of the resistances due to the two systems, ignoring any 
interference effects. 



FIG. 2 




The graph is of interest in its exhibition of the typical humps and 
hollows, occurring in general when 27rc^/^ is a sub-multiple of 2h. The 
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prominence of the interference effects depends upon the relative magnitudes 
of the constants / and A ; the example we have chosen shows a pronounced 
effect due to the final maximum of interference being near the maximum of 
the mean curve. 

8. We may note briefly the interpretation in terms of a submerged body 
The surface of the body is one formed of stream-lines due to the two equal 
doublets in a uniform stream ; the axes of the doublets are in the same 
horizontal line at a depth /, which must be large compared with, at least, 
the vertical dimensions of the body. For instance, with suitable relations 
between the constants, the result would give the wave resistance of two 
small bodies, of nearly spherical shape, one behind the other at a distance 
large compared with their dimensions. 

9. By combining simple symmetrical pressure systems, we may generalise 
the previous results; this seems an easier process than the direct discussion 
of unsymmetrical systems. We shall assame that the component simple 
systems are all of the type (12) and have the constant / of the same value, 
and that the centres of the systems all lie on the axis Ox, 

In the first place we must extend the analysis of §8 to two components of 
unequal magnitudes A and B, with their centres at the points (A, 0) and 
(Jc, 0) respectively. From the argument expressed in (29)--(32), it is easily 
shown that the value (33) for the wave resistance must be altered by replacing 
A^ cos^ (fcoh sec (j>) by 

i[A2 + B2 + 2ABcos {/^o(A-^)seo^}]. (34) 

Suppose further that the pressure system is given as a continuous line 
, distribution of components along Ox in a range from hi to As, the magnitude 
of the element with its centre at (x, 0) being proportional to some function 
-v/f (x) ; in other words, suppose the surface pressure is given by 



j? = A/ 



rh 



-ifr (h) dh 



the function '\^ (A) being such that the transformations used in the preceding 
analysis are permissible. For the system (35) we have to sum (34) for all 
possible pairs of elements ; this is performed by taking the double summation 



1 



4 



h. 



-^ Qh) dh '^ (k) cos {fco Qi-^h) sec ^} dh, (36) 



The wave resistance for the system (32) can be completed now from (33) 
and (36) ; we have 



'hj, thn 



f (A) dh ^-f (h) dh see/'^ ^e-2Ko/sec2.|> ^os {kq (h-^h) sec ^} dcj). (37) 

VOL. XOV.— A. 2 G 
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10. Consider as an example the case when yfr{h) is constant, so that the 
surface pressure is 



p/Af = f - 






^ This may be regarded as the combination of two equal systems of opposite 
sign, with their centres at the points (h, 0) and {--Ih 0) but not symmetrical 
round these points. 

In this case, after carrying out the integrations with respect to h and k 
(37) gives 

R = (16 nr/gp) A^/co | "" sec^ (j^e-^'ofBecH gin^ (fcoh sec <^) dcj). (39) 

Jo 

The integral may be treated similarly to (33). One of the main differences 
lies in the factor sin^ (/^o^secc^) instead of cos^ (a:o^?' sec ^) ; this is because 
we have now two equal positive and negative systems instead of two positive 
systems, and in consequence the series of humps and hollows on the resistance 
curve will be interchanged. 

We have chosen this case partly because of the corresponding problem in 
the motion of a submerged body at depth /. Integrating a line of doublets 
of constant strength results in a simple source at one end of the line and an 
equal sink at the other. Hence, the submerged body is one of the oval- 
ishaped surfaces of revolution formed by combining a source and sink with a 
uniform stream ; it follows that, as in § 5, the strength of the source is 
A.GJ2gp, It may be noted that the coefficient A in (39) has different 
dimensions from that in (33), agreeing with its introduction in (38). By 
making h small in (39) we regain the former result for a sphere. 

11. If a prolate spheroid of semi-axis a and eccentricity e is moving in an 
infinite liquid with velocity c in the direction of its axis of symmetry, it can 
be shown that the velocity potential may be written in the form 

^ "" J -ae {(x^hf + f + Z^Yl'' ^^^^ 

where A = l/[46/(l— e^)— 21og {(l + e)/(l— ^)}1 and where we have, for the 
moment, taken Ox along the axis of symmetry of the spheroid. This 
expresses ^ as due to a line of doublets ranged along the axis between the two 
foci. Hence the surface pressure corresponding to the motion of the spheroid 
with its axis at depth/ is 

7,\2 I ...2 I .2^ S/2' V^-*-/ 



p - 2^p/A j_^^^^^^^^_^^^,^^,^^,^3/„ 
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reverting to axes with the origin in the free surface. It should be noticed 
that, as in § 5, the surface pressure (41) does not give the same surface 
elevation as the moving spheroid; the surface condition in the latter case is 
that the pressure should be constant at the free surface. But (41) does give 
the same wave formation as the spheroid, and that is the part of the surface 
effect upon which the wave resistance depends. The complete surface eleva- 
tion can be easily written down by direct methods as in the case of the 
submerged sphere. 

Using (41) now as an example of (35), we find the wave resistance of the 
spheroid from (37) ; after integrating with respect to h and h, the result is 

E = USiT^gpah^A^ r'sec2^e-2«o/«ec^<^{j3^2(/co^^sec^)}2<^^. (42) 

It can be verified that this gives the result for the sphere by making e zero. 

Eor a given relation between/ and ae, the value of E can be obtained 
approximately by the numerical methods used in the previous examples ; 
judging from rough calculations, it appears that the resistance curve does not 
show prominent humps and hollows. This might be anticipated from the 
surface pressure (41), which can be evaluated in simple form ; if we repre- 
sent the pressure distribution by a surface with p, x, and y as co-ordinates, 
then (41) gives a single oval-shaped peak with its longer axis in the direction 
Ox. On the other hand, the pressure distribution (28) represents two 
distinct peaks. We may compare in this respect the behaviour of ships' 
models ; it depends upon the shape of bow and stern, and the relation 
between them, whether the resistance curve has marked interference effects 
or is a continuously ascending curve. 

12. We have limited the previous cases to combinations of simple pressure 
systems ranged along the axis Ox. The method can obviously be extended 
to systems with their centres on Oy ; or again, for systems situated in the 
plane ^3/^, a four-fold summation in the manner of (36) would give further 
generality. For the corresponding problem of the motion of a submerged 
body, one could obtain the wave resistance of any body whose surface is 
formed of stream lines due to the combustion of sources and sinks with a 
uniform stream. 
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